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Quenching across quantum critical points: role of topological pat- 
terns 
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Abstract. - We introduce a one-dimensional version of the Kitaev model consisting of spins on 
a two-legged ladder and characterized by Z2 invariants on the plaquettes of the ladder. We map 
the model to a fermionic system and identify the topological sectors associated with different Z2 
patterns in terms of fermion occupation numbers. Within these different sectors, we investigate 
the effect of a linear quench across a quantum critical point. We study the dominant behavior of 
the system by employing a Landau-Zener-type analysis of the effective Hamiltonian in the low- 
energy subspace for which the effective quenching can sometimes be non-linear. We show that 
the quenching leads to a residual energy which scales as a power of the quenching rate, and that 
the power depends on the topological sectors and their symmetry properties in a non-trivial way. 
This behavior is consistent with the general theory of quantum quenching, but with the correlation 
length exponent v being different in different sectors. 



Of late, two different concepts have instigated a surge 
of active research in quantum many-body phenomena - 
the notion of topological order p}j4] and quenching across 
quantum critical points (QCPs) [MS]. The former, de- 
scribed by global invariants, has been keenly studied for 
its fundamental significance, realizations in physical sys- 
tems, associated phase transitions and potential applica- 
tions to quantum computation. The latter has offered a 
mine of valuable information on the nature of the QCP 
in question, for instance, in the scaling behaviors of the 
defect density and the residual energy upon quenching. 
These quantities are governed by a quantum version of 
the Kibble-Zurek mechanism and in d dimensions, they 
scale with the quenching rate 1/t as 1 / T <WO iy + 1 ) ; where 
v and z are respectively the correlation length and dynam- 
ical critical exponents of the QCP. Post-quench dynamics 
is thus effectual in extracting a combination of these ex- 
ponents defined by the growth of the correlation length 
as I A — A c | _l/ and the relaxation time as |A — X c \~ zv , as a 
parameter A in the Hamiltonian of the system approaches 
a critical value A c 0[6]. In this Letter, we demonstrate 
that a marriage between these two concepts makes for a 
synergistic union. Specifically, we show that invariants in 
a system having topological order can impose dramatic 
constraints on statics and dynamics, giving rise to quali- 



tatively different post-quench behaviors; hence, quenching 
can act as a powerful means of discerning patterns and 
symmetries in the topological invariants of a system. 

Our focus here is on models sharing traits of the cele- 
brated Kitaev model [3], namely spin- 1/2 lattice systems 
with an exponentially large number of configurations dis- 
tinguished by conserved Z2 quantum numbers, which we 
refer to as topological sectors. Detailed studies charac- 
terizing the various topological sectors and distinguishing 
their physical effects are sparse; here we address these is- 
sues in the context of quenching in a simple two-legged 
ladder model which captures salient features of the two- 
dimensional parent Kitaev system. We explicitly char- 
acterize the Z 2 patterns based on their periodicity and, 
within these different sectors, study quenching at zero 
temperature through a QCP. We find that the residual 
energy, the difference between the post-quench energy ex- 
pectation value and the true ground state energy, respects 
the dv j(zv + 1) power-law scaling form with respect to 
the quench rate. While the dynamical exponent takes the 
value z = 1, the value of v, and therefore the power-law, 
depends on the topological sector being probed. For sec- 
tors endowed with higher symmetry in their periodicity, 
v deviates from unity and the power-law is no longer of 
the 1/t 1 / 2 form typical of many one-dimensional systems 
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n-l,b Ji n,b n+l,b n+2,b 

Fig. 1: Picture of the two-legged Kitaev ladder showing the 
couplings Ji and the ^-valued operators W n - 

[5,6 . We thus theoretically demonstrate, using a simple 
model, the effect of topology on quantum critical behavior 
and the rich resultant physics in quench dynamics. 

The quasi-one-dimcnsional Kitaev model that we con- 
sider here (see Fig. [IJ consists of a two-legged ladder with 
a spin-1/2 at each site. We label the sites as (n, a) and 
(n, b), where a, b denote the upper and lower legs respec- 
tively. The Hamiltonian for the system is given by 

H = £[-A««Al,a + <6^+l,&) + ^<,a<,bl W 
n 

where <r" a ^ b denote the Pauli matrices; the couplings con- 
sidered here are analogous to those associated with a strip 
of the honeycomb lattice in the parent Kitaev model. (We 
set h = 1.) We assume J\ > 0; negative Ji values can 
be gauged away by performing 7r-rotations of appropriate 
spin components. As discussed in Ref. [10], the system 
exhibits a rich phase diagram as a function of the cou- 
plings Ji and J3. Each plaqucttc of the ladder hosts a 
topological invariant and together these invariants char- 
acterize different topological sectors. For a given pla- 
quette (see Fig. [1]), we can define a Hermitian operator 

W n = ' 7 n,a cr ^6 cr n+l,a cr ,1+l,6- 0ne Can SnOW tnat the W n 'S 

commute with each other and with the Hamiltonian, and 
that W% = 1. The invariant eigenvalues, W n = ±1, pro- 
vide a set of Zi quantum numbers and, for a system with 
N plaquettes, yield 2 N distinct sectors corresponding to 
different sets of these numbers. 

As in previous studies [T0J[TT], we map the system 
to that of two sets of spinless fermions, one probing 
the dynamics and the other associated with the Kitaev 
model's hallmark Ising anyons pQ. The mapping in- 
volves a Jordan- Wigner transformation which takes the 
form a n = S n , a c n = S n ^ a cj? Q , b n = S n> h b , 

d n = S nib a v n b , where S n>a / h denote a string of cr na/b 's 
defined as follows. Assuming that each leg contains N 
sites with open (rather than periodic) boundary condi- 
tions, S rit a = llm=i a m a weaves from the left to the right 
along the top leg in Fig. [T]to the (n— l)-th site, and S n .b = 
°m b weaves completely along the 
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top leg left to right and then along the bottom leg from 
right to left to the (n + l)-th site. The transformed Hermi- 
tian operators are Majorana fermion operators satisfying 
relations such as a n = 1, {a m ,a n } = 25 mn , {a m ,b n } = 0, 
and so on. The resultant Hamiltonian takes the form H = 



E„[-^i( » a »+i + Kb n +i) + j(-l)"J3a n () n n* =n W m ], 
where W n = —{ic n d n )(ic n +id n +i). Thus, the Hamilto- 
nian is quadratic in the a — b fermions and dependent on 
the Z2 values of the topological invariants. 

Next, the Majorana fermions can be combined pairwisc 
to form Dirac fermions, providing an understanding of 
the system in terms of Dirac fermion occupation num- 
bers. The natural combinations, which pair partners on 
each rung of the ladder, are /„ = (i n /2)(a n + ib n ) and 
g n = (l/2)(c„ + id n ); these satisfy {f m , f£} = {g m ,gX\ = 
8 m ^ n - The topological invariants can be expressed in terms 
of occupation numbers of the g-fermions {g\g n = 0, 1) as 
W n = —s n s n+ i where s n = ic n d n = 2g} l g n - 1. The 
value of W on the n-th plaquette is +1 if its bordering 
g- fermions have opposite occupation numbers and —1 if 
they are the same. The s n 's themselves are non-local in 
terms of the W„'s: s„ = Il^n^ 7 ™- Specifying the val- 
ues s n — ±1 completely fixes those of the W n 's. For in- 
stance, alternating values of s n along the ladder gives all 
W n = +1, a constant value of s n gives all W n — — 1, and a 
domain wall in s n is equivalent to embedding a W n = +1 
at the wall amidst all other W n 's being — 1. [Changing 
the values of W n and thus the topological sector can be 
achieved by the appropriate local action of operators such 
as (7% a cr„ +1 a and a v n b (Jn+\ bi which can alter the occupa- 
tion numbers of the g-fermions on any given site]. In the 
thermodynamic limit, the Hamiltonian in Eq. ([IJ takes 
the local form 

H = £>2Ji(4 +1 /„ + flf n+1 ) + J 3 s„(2/t/„ - 1)], (2) 

n 

namely a tight-binding system for the /-fermions whose 
on-site chemical potential is determined by the g-fcrmion 
occupation number. The Hilbert space of 2 2N associated 
with the 2N spins on a 7V-plaquette ladder is thus split 
into two Hilbert spaces, each of dimension 2 N , correspond- 
ing to those of the 'dual' / and g fermions (whose roles can 
be exchanged by swapping the x and y spin components.) 
The /-fermions probe the dynamics of the system, while 
the g-fcrmions, whose Hilbert space dimension reflects the 
SU (2)2 structure associated with Ising anyons pQ, deter- 
mine its topology. 

The topological sectors for this model can be identified 
by the patterns exhibited by the s„'s. For patterns of s n 
having a period p, the Hamiltonian of Eq. ([2J splits up 
into decoupled sub-systems involving p momentum states 
with values k + 2-nq/p, where q is an integer satisfying 
< q < p — 1 and k lies in the range — ir < k < —tt + 
2n/p. To determine the filling of the /-fermions for the 
lowest energy state within a given sector, consider the limit 
J 1 = and J3 < 0. Here energy minimization requires 
that the occupation number f\fn at site n be equal to 
(1 + s n )/2; the gap to the first excited state is given by 
2|J 3 |. A periodicity having m of the s„'s be +1 and p — m 
be —1 thus has m sites occupied and the others unoccupied 
within a period; the filling is m/p. Turning on J\ does 
not change the particle number; hence the filling remains 
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constant in the lowest energy state. As can be shown 
perturbatively around the J\ = limit, a small value of 
Ji changes the ground state energy, to second-order in J\ , 
by AE = -(J1/J3) Ytni 1 ~ s nSn+i)- Hence, the ground 
state sector in which the energy of the lowest-lying state 
is the minimum amidst all the sectors is the one in which 
s n = (-1)" and W n = 1 for all n. 

Having characterized the system by the different topo- 
logical sectors and their eigenstate structures, we now in- 
vestigate the effects of quenching within different topolog- 
ical sectors. We assume that the system is initially placed 
in a particular sector; a key feature of a topological sector 
is that the system will remain in that sector at all times if 
there are no perturbations which can dynamically change 
the topological invariants. We initialize the system in its 
lowest energy state at J3 = — 00; at that point, the en- 
ergy gap between the ground state and the excited states 
is infinitely large. We then vary J3 linearly in time as 
Js(t) — Jit/r; this quench tunes the system through a 
QCP at J 3 = as discussed below. We quantify the effect 
of the quench via the residual energy per site attained at 
the final time, which is defined as 



E r — 



lim 

t,JV->-o 



(H) f - E 

N\Mt)\ ' 



(3) 



where (H) / denotes the expectation value of the Hamilto- 
nian in the final state reached, and Eq is the ground state 
energy at t — > 00. We study the dependence of E r on r for 
J\T 3> 1; this must vanish in the adiabatic limit r — > 00. 
We find that E r shows different qualitative trends depend- 
ing on the manner in which the topologically sensitive J3 
term in Eq. ([2]) couples different momentum modes. We 
discuss three distinct classes of behavior: (i) the trivial 
case of decoupled modes, (ii) a Landau-Zener type r -1 / 2 
dependence of E r due to direct coupling between the low- 
energy modes, and (iii) the most interesting case of a non- 
trivial power-law due to an indirect coupling. 

Let us represent the /-fermions in the momentum basis 
as f k = ^Ejne _!t ", where -tt < k < tt. In this 
basis, the sector having all W n = — 1, i.e., all s n = 1 for 
all n, can be seen to behave trivially under quenching. The 
modes having different momenta do not mix for any J3 in 
the Hamiltonian of Eq. Thus, if one begins in the 

lowest energy state for J3 = —00, i.e., with fXfn = 1, one 
remains in that state as J3 changes to 00 where it becomes 
the highest excited state. Hence, the residual energy E r 
is equal to 2 independently of the value of r. 

The ground state sector having all W n = 1 offers an in- 
stance of direct coupling between low-energy modes. As- 
suming that s n = (—1)", the Hamiltonian decouples into 
sub-systems having pairs of momenta k and k + tt as 
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4 Ji cos k 
2J 3 



2J 3 
-4 Ji cos k 



k2- 



(4) 



where fk2 = (fk+-w,fk) T represents the momentum mode 
annihilation operators (T denotes the transpose), and 



— tt < k < 0. The corresponding eigenenergies are 
E k ± = ±-\/4J| + 16 J 2 cos 2 k. For J 3 = 0, the energies 
vanish at k — — ir/2, indicative of a quantum critical point 
(QCP); in fact, a QCP having gapless modes at k c = ±tt/2 
occurs at J3 = in any topological sector with a half- 
filled ground state. Here, the dynamical critical exponent 
is z = 1 and the correlation length exponent is v — 1 since 
the energy vanishes as \k — k c \ at the QCP at J3 = and 
as I J3I at k = k c - 

For the quench J3 = Jit/r, a 7r/2 unitary rotation inter- 
changes the diagonal and off-diagonal terms in the matrix 
appearing in Eq. (fj|, exactly mapping the quench to the 
well-known Landau-Zener problem j6l[7l[T2]. The proba- 
bility pk of ending in an excited state at t — > 00 and the 
net residual energy, for which each sub-system contributes 
4pk, are therefore given by 



Pk 



exp[— 8ittJi cos 2 k], E r = 



dk 
2^ 



4p fc . (5) 



The probability pk is largest for the low-energy modes near 
k c = — 7r/2 since this is where the gap between the two 
states vanishes for J3 = 0. In the limit tJ\ — > 00, the 
residual energy is dominated by this low-energy regime; 
expanding around k c yields a Gaussian integral and the 
power-law form E r ~ 1/r 1 / 2 . This scaling is consistent 
with the values of the critical exponents given above and 
in keeping with quenches through QCPs in many one- 
dimensional systems [5]E|- 

We now turn to an instance of a topological sector yield- 
ing a completely different quench power-law of 1/r 2 / 3 
due to a higher symmetry in its periodic structure, as 
we show by heuristic arguments and numerics. The in- 
stance is of W n — (—1)™, or equivalently, the signs of s n 

forming the half-filling period 4 pattern (+ H ). The 

pattern can be expressed as s„ = — \/2 cos(7m/2 + 7r/4). 
The resulting Hamiltonian of Eq. ^ decouples into sub- 
systems involving four momenta with annihilation opera- 
tors f k 4 = (fk+3n/2>fk+w, fk+ir/2,fk) T , where k lies in the 
range -tt < k < -tt/2. Explicitly, H 4 = Y,k Jki H kifk4, 
where Hki = Mk + Nk with 



M k = -J3V2 
N k = 4Jr 



f3 
(3* 

— sin k 

cos A: 



f3* 
13 



(7) 



Here (3 = e l7r / 4 , /i Q denote Pauli matrices and / is the 
identity matrix. The higher symmetry of this sector is re- 
flected by the fact that while the four momentum states 
cannot be decoupled into pairs, they are not all directly 
coupled to one another. The eigenenergies of HkA come in 
pairs ±Ek4, as can be derived from the symmetry property 
UHkiU^ = -HkA, where U = <g> p x . The QCP asso- 
ciated with this system occurs, as reflected in Eqs. (J6JTJ) , 
at J3 = and k close to —tt and —tt/2. For a quench 
J 3 = Jit/r, the dominant contribution to the residual 
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energy comes from particle-hole pairs closest to the Fermi 
level at zero energy for which the cost of quenching into an 
excited state is low. For instance, for momenta k gl — tt, 
the relevant pairs closest to the Fermi energy at small 
J3 are |1,0) and |0,1), where n p and no in state \n p ,no) 
denote /-fermion occupation numbers of momenta k and 
k+w, respectively. These states are not directly coupled by 
the Hamiltonian Hk4 and only mix via the higher energy 
states at k + tt/2 and k + 3n/2. To second-order pertur- 
bation in J 3 , this coupling is of order Jf/Ji; the effective 
Hamiltonian for the two-level system, in the appropriate 
basis, takes the form 

Hk4,eff = Jl ( 1 ll _f jT 2 ) , (8) 

where we have expanded — sin(fc — tt) ~ fc, and have used 
the time-dependent form of J3 = J\t/r. A similar form 
applies to modes at k S —tt/2. Arguments similar to the 
previous case show that here z = 1 and v — 2. 

The effective quench in Eq. © is quadratic in time. 
Unlike the linear quenching problem, no exact solution is 
known for the excitation probability pk for the Schrodinger 
equation idipk/dt = Hk4, e ffil>k given by Eq. ©. However, 
we can invoke scaling arguments to find the power-law de- 
pendence of pk on r. Rcscaling time as t' = t/r 2 ^ 3 results 
in the excitation probability being governed by the pa- 
rameter fcr 2 / 3 . If fcr 2 / 3 < 1, the corresponding modes 
have a significant weight for occupying the excited state 
as t' — > 00, while the modes with fcr 2 / 3 3> 1 remain in 
the ground state. Integrating pj~ around the low-energy 
modes to obtain the residual energy thus results in the 
scaling E r ~ 1/r 2 / 3 . We ascertain this behavior by nu- 
merically solving the time-dependent Schrodinger equa- 
tion given by the full-fledged Hamiltonian in Eqs. ([B][7]) 
from t = —00 to 00 with J3 = J\t/r. We start in the 
lowest energy state at J3 = —00 consisting of the two 
occupied states (1, — 1, 1, — 1)/2 and (1, i, — 1, — i)/2, and 
time evolve to obtain the probability pk for occupying the 
excited states. The left panel of Fig. [5]shows the resultant 
Pk versus the scaled variable (fc + 7t)(Jit) 2 / 3 for k — tt; 
the curve is independent of r for J\T 3> 1. The probabil- 
ity pk near k ;$ —tt/2 is related to this curve by mirror 
symmetry. Analogous to Eq. (|5|) , integrating 2pk /tt over k 
from — tt to —tt/2 yields the residual energy E r . The right 
panel of Fig. [2] shows a plot of log(-E r ) versus log(Jir). 
The clean linear fit with a slope close to —2/3 confirms 
our predicted 1/r 2 / 3 scaling form. 

This unusual power-law scaling is generic to a series of 
topological sectors that are at half-filling, i.e., those that 
in a period 2m have half the s„'s take on the value +1. 
Of the associated 2m momentum modes, the states clos- 
est to k = ±tt/2 and thus to the Fermi energy dominate 
the quench. As in the example above, the quench through 
the QCP at J 3 = shows unusual scaling if these states 
are not directly coupled. This corresponds to the J3 term 
in the Hamiltonian having no matrix element between the 
low-energy modes or, equivalcntly, having A m vanish in 




Fig. 2: Left panel: Plot ofpj; versus (k+n)( Jit) 2 ^ 3 for k ^ — n, 
and Jit = 2 (red +), 8 (black x) and 32 (blue *). Right 
panel: Logarithmic plot of E r versus J\t in the sector with 
W n = (-1)". A linear fit gives E r = 0.14/( Jir) ' 67 which is 
close to a —2/3 power-law. 

the Fourier expansion At = Y^j=i Sje l7rj7 / m . Clearly the 
number of such cases increases with increasing m. For 
instance, such examples of period 8 which also yield the 
1/t 2 / 3 power-law by way of low-energy modes being cou- 
pled via one intermediate high energy state are the sets 

(+ + + H ), (+ + + 1 ), (+H h + ), 

(+H 1 1 — ), and cyclic permutations thereof. How- 
ever, in cases of even higher symmetry where the low- 
energy modes are only coupled via a path involving q inter- 
mediate high energy states, the effective coupling is of the 
form Ji(t/r) q+1 for the linear quench J3 = Jit/r. In prin- 
ciple, the generalized scaling argument presented above 
would predict that the residual energy would have the 
scaling form E r ~ T -(i+ 1 )/(i+ 2 ) corresponding to z = 1 
and v = q + 1, but the analysis would depend on the spe- 
cific topological sector and dominant couplings close to 
the QCP. 

In conclusion, we have demonstrated that in certain sys- 
tems, topological order can play a major role in deter- 
mining the universality class and behavior close to QCPs. 
Topology can greatly constrain dynamics and yield un- 
usual power-law scaling in quenches across QCPs which 
can distinguish different orderings based on symmetry 
properties. The one-dimensional example showcased here 
ought to be realizable using recent schemes proposed in 
cold atomic systems [T3] and the associated topological 
sectors ought to be accessed using local operations. We 
expect that topology-driven scaling forms would be man- 
ifest in several observables such as dynamic spin-spin cor- 
relations, and in two-dimensional generalizations of the 
Kitacv model. Finally, due to a finite quantum critical 
regime in parameter space and the robustness of topolog- 
ical sectors against a large class of perturbations, we ex- 
pect the features discussed here to persist at temperatures 
much smaller than the energy gap at the initial time. At 
higher temperatures, the interplay between thermal and 
quantum fluctuations on quench dynamics |14| would re- 
quire further investigation in our model. 
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